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Abstract. We present a new study of the high latitude galactic contributions to the millimeter sky, based on an 
analysis of the WMAP data combined with several templates of dust emission (DIRBE/COBE and FIRAS/COBE) 
and gas tracers (HI and H Q ). Fo study the IR to millimeter properties of the diffuse sky at high galactic latitude, 
we concentrate on the emission correlated with the HI gas. We compute the emission spectrum of the dust/free- 
free/synchrotron components associated with HI gas from low to large column densities. A significant residual 
WMAP emission over the free-free, synchrotron and the dust contributions is found from 3.2 to 9.1 mm. We 
show that this residual WMAP emission (normalised to 10 20 atoms/cm 2 ) (1) exhibits a constant spectrum from 
3.2 to 9.1 mm and (2) significantly decreases in amplitude when Nhi increases, contrary to the HI- normalised 
far-infrared emission which stays rather constant. It is thus very likely that the residual WMAP emission is not 
associated with the Farge Grain dust component. Fhe decrease in amplitude with increasing opacity ressembles 
in fact to the decrease of the transiently heated dust grain emission observed in dense interstellar clouds. Fhis 
is supported by an observed decrease of the Hi-normalised 60 fim emission with HI column densities. Although 
this result should be interpreted with care due to zodiacal residual contamination at 60 /im, it suggests that the 
WMAP excess emission is associated with the small transiently heated dust particles. On the possible models of 
this so-called "anomalous microwave emission" linked to the small dust particles are the spinning dust and the 
excess millimeter emission of the small grains, due to the cold temperatures they can reach between two successive 
impacts with photons. 
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1. Introduction 

At millimeter wavelengths, one of the major challenges in 
high sensitivity Cosmic Microwave Background (CMB) 
anisotropy study is to determine the fraction of the 
observed signal due to diffuse galactic foregrounds. 
Three different components have been firmly identified 
at high latitudes >10°): thermal dust emission, 

synchrotron and free-free. Dust emission dominates the 
far-infrared surveys. Its spatial distribution and frequency 
dependence are quite well-determined for wavelengths 
shorter than ~800 /im. Above ^800 /im, present data 
currently do not give any strong constraints. So far, dust 
emission estimates in the millimeter range are thus an 
extrapolation of what is known at shorter wavelengths. 
Synchrotron radiation dominates radio-frequency sur- 
veys, but Banday & Wolfendalc (1991) and Bennett et 
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al. 1|1992[) showed that the spectral index steepens with 
frequency and exhibits spatial variations which are poorly 
known. Free-free emission has a well-determined spectral 
behavior and templates are now available thanks to the 
WHAM H Q survey of the northern sky (Reynolds et al. 
rreggl Haffner ITMty and the SHA SSA H Q survey of the 
southern sky (Gaustad et al. 12001(1 . 

Cross-correlations of CMB data with far-infrared 
maps have revealed the existence of a microwave emis- 
sion component (the so-called "anomalous microwave 
emission") with spatial distribution traced by these 
maps. This component has a spectral index suggestive 
of free-free emission and so has been first interpreted as 
free- free emission (Kogut et al. 119961 . However, Kogut 
( 1999) showed in small parts of the sky that were covered 
by H a data that the microwave emission was consistently 
brighter than the free- free emission traced by H Q . Thus, 
the correlated component cannot be due to free-free 
emission alone. This is confirmed more recently by 
Banday et al. PUU3|) also using COBE/DMR data. 

Recent works suggest that this anomalous far-infrared 
correlated component originates from spinning dust grain 
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emission (Draine & Lazarian 1998a, Dc Oliveira-Costa et 
al. HMfl Effing , tentatively detected at 5, 8 and 10 GHz 
by Finkbeiner et al. (2()02i. An alternative explanation 
is provided by thermal fluctuations in the magnetization 
of interstellar grains causing magnetic dipole radiation 
(Draine & Lazarian 1999). However, very recently, 
Bennett et al. (|2003|) using WMAP data do not find any 
evidence for the anomalous microwave emission. Their 
foreground component model comprises only free-free, 
synchrotron and thermal dust emission, and the observed 
galactic emission matches the model to <\%. Note that in 
their global analysis, they are dominated by the brightest 
parts of the sky i.e. the galactic plane and the high 
latitude dense interstellar clouds. Thus, results may not 
apply to the most diffuse regions. 

We present in this paper a new study of the galactic 
contributions to the millimeter sky, based on an analy- 
sis of the WMAP data combined with several templates 
of dust emission (DIRBE/COBE and FIRAS/COBE) and 
gas tracers (HI and H Q ). We focus only on the high lat- 
itude regions where the results are easier to interpret in 
term of physical properties of dust and where CMB anal- 
ysis are performed. The paper is organised as follows. We 
first present the data we use together with their prepara- 
tion (Sect.0). We then derive the spectrum (from 100 /jm 
to 10 mm) of the Hl-correlated component (Sect 3.1) 
and show that there exists a residual microwave emission 
(over free-free, synchrotron and far-infrared dust emis- 
sion) whose Hl-normalised amplitude decreases when the 
HI column density increases but without any significant 
spectral variations (Sect. 3.2). We then discuss the results 
in Sect. 4. 

2. Data-sets used and their preparation 

2.1. CO BE data 

The COBE satellite was developed by NASA's Goddard 
Space Flight Center to measure the diffuse infrared and 
microwave radiation from the early universe to the lim- 
its set by our astrophysical environment. It was launched 
November 18, 1989 and carried three instruments: 

— a Far Infrared Absolute Spectrophotometer (FIRAS) 
to compare the spectrum of the cosmic microwave 
background radiation with a precise blackbody (at 
each sky position, with an angular resolution of 7°, 
we have one spectrum from 1 to 97 cm -1 ) 

— a Differential Microwave Radiometer (DMR) to map 
the cosmic radiation (at 31, 53, 90 GHz with a 7°beam) 

— a Diffuse Infrared Background Experiment (DIRBE) 
to search for the cosmic infrared background radiation 
(10 photometric bands from 1 to 240 /im with an an- 
gular resolution of 40'). 

COBE data are presented in a quadrilatcralized spher- 
ical projection (the so-called COBE Quadrilateralized 



Spherical Cube, CSC), an approximately equal-area pro- 
jection in which the celestial sphere is projected onto an 
inscribed cube. The DIRBE convention is to divide each 
cube face into 256x256 pixels; thus all sky-maps have 
256x256x6 = 393216 pixels. Each pixel is approxima- 
tively 0.32°on a side. For FIRAS and DMR, each cube face 
is 32x32 pixels leading to a total of 6144 pixels (of ~2.6°). 

We use the so-called (i) "Sky Maps and Analyzed 
Science Data Sets" DMR Data (ii) "Galactic Dust 
Continuum Spectra and Interstellar Dust Parameters" 
FIRAS data, that give the residual sky spectrum after 
modelled emission from the CMB, zodiacal emission, 
and interstellar lines have been subtracted, (hi) DIRBE 
"Zodi-Subtracted Mission Average (ZSMA) Maps" for 
which the zodiacal light intensities were subtracted 
week by week and the residual intensity values were 
averaged to create Maps. All COBE data are availabe at 
http: / /lambda. gsfc.nasa.gov/product / cobe 



2.2. Gas tracers 

The HI data we used are those of the Leiden/Dwingeloo 
survey which covers the entire sky down to S = —30° with 
a grid spacing of 30' in both 1 and b. The 36' half power 
beam width of the Dwingeloo 25m telescope provides 21 
cm maps at an angular resolution which closely matches 
that of the DIRBE maps. Details of the observations and 
correction procedures are given by Hartmann (1994) and 
by Hartmann & Burton ( 1997). It should be noted that in 
this data-set special care was taken for the removal of far 
sidelobes emission which makes it particularly suitable for 
high latitude studies. We derive the HI column densities 
using 1 K km s -1 =1.82 10 18 H cm~ 2 (optically thin 
emission). 

Thanks to the WHAM survey of the northern sky 
(Reynolds et al. IT7M1 Haffner ITMfll and t he SH ASSA 
survey of the southern sky (Gaustad et al. I2()()l|l . it is 
now possible to have a whole sky map of the H Q emis- 
sion (Dickinson et al. 1277051 Finkbeiner l2UU3|) . Since the HI 
maps cover the sky down to S = —30°, the H a emission 
we use is our analysis is mostly given by the WHAM sur- 
vey. WHAM provides a 12 km s" 1 velocity resolution with 
one-degree angular resolution down to sensitivity limits of 
0.2 R (1 R = 10 6 /4?r ph cnT 2 s" 1 sr" 1 ) in a 30 second ex- 
posure. The one-degree angular resolution nicely matches 
the DIRBE resolution. We use the H Q map and the conver- 
sion factors to free-free emission (using T c =8000 K) from 
Finkbeiner (2003) to derive templates of free- free emis- 
sion. Since we work only on high latitude regions, the H Q 
emission has not been corrected for extinction (the dust 
absorption is likely to be very small, less than 5%). The 
free-free templates are used to derive a well-understood 
contribution to the millimeter channels. 
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2.3. Synchrotron templates 

Synchrotron emission arises from relativistic cosmic ray 
electrons spiralling in the galactic magnetic field. This 
emission dominates surveys at radio frequencies. The only 
all-sky map at low frequencies that probe the synchrotron 
emission is the 408 MHz survey of Haslam et al. (1982). 
For many years, this map has been used to predict the 
synchrotron emission in the millimeter channels, assum- 
ing a frequency dependence with a constant spectral in- 
dex of about 2.75. However, Bennett et al. 1|2003[) have 
shown that the synchrotron spectral index exhibit strong 
spatial variations and is steeper in the WMAP bands than 
at radio frequencies. We thus use the WMAP synchrotron 
maps derived by Bennett et al. (J2003 ) using the Maximum 
Entropy Method as frequency dependent well-determined 
synchrotron templates. 

2.4. WMAP data 

The WMAP 1 (Wilkinson Microwave Anisotropy Probe) 
mission, designed to determine the geometry, content, and 
evolution of the universe, has successfully provided full sky 
maps at 23, 33, 41, 61 and 94 GHz at respectively 13.2, 
39.6, 30.6, 21 and 13.2' FWHM resolution with unprece- 
dent sensitivity. A detailed description of the delivered 
data-sets for the first 12 months of operation of WMAP is 
given in the WMAP: explanatory Supplement (editor M. 
Limon et al., Greenbelt, MD: NASA/GSFC). The data 
we use are the first-year "Smoothed I maps" which are 
the temperature maps at each frequency, smoothed to a 
common resolution of 1 degree. Data are delivered in the 
HEALPix 2 format with N sidc =256 



is locally heated by nearby stars (like the HII regions) fol- 
lowing Lagache et al. l|1998f) . We stress out that this latter 
pixel selection, although necessary to keep in the analysis 
only diffuse parts of the sky, does not change the results 
and conclusions of the paper. 



3. Analysing the data 

3.1. Deriving the Hl-correlated component spectrum 

The Hl-correlated dust emission is the dominant 
component at high galactic latitude at infrared/far- 
infrared/submillimeter wavelengths (except in the very 
low HI column density regions where the Cosmic Infrared 
Background becomes an important contribution, e.g. 
Lagache et al. I1999fl . One way to study the infrared to 
millimeter properties of the diffuse sky at high galactic 
latitude is therefore to concentrate on the emission 
correlated with the HI gas. We search here for the 
spatial/spectral variations of the infrared to millimeter 
emission with the HI gas column densities. 

To compute the emission spectrum of the component 
associated with HI gas from low to large column densi- 
ties, we use a differential method that removes, within 
statistical variance, any residual infrared emission that is 
not correlated with the HI gas such as an isotropic com- 
ponent. We first select sky pixels according to their HI 
column density and sort them into sets of pixels brack- 
eted by selected values of Nhi- Correlated HI emission 
spectra are then computed for each set of pixel k using 
the equation: 



2.5. Data preparation 

All the data have to be put in the same projection and 
at the same resolution. The resolution is set by the 
FIRAS experiment since it is the lowest resolution of our 
data-sets (7°). All data but DMR and FIRAS are thus 
converted in the DIRBE CSC format and then convolved 
with the FIRAS beam and degraded to the FIRAS CSC 
resolution (see Lagache 120031 for more details). 

We have removed for each data-set the cosecant law 
variation (1) to avoid the obvious large scale correlations 
between all galactic components concentrated in the 
disc and (2) to be consistent with the WMAP data that 
measure only differentially on the sky and thus does not 
measure the largest angular scales. 

We restrict our analysis to \b\ > 15° and exclude the 
Small and Large Magellanic clouds, together with the p - 
Ophiucus complex. We also remove cold molecular com- 
plexes (as the Taurus cloud), and regions where the dust 



F v (k) = 



< F > fc - < F >p 
< Nhi >k - < N HI > 



(1) 



1 http:/ /lambda. gsfc.nasa.gov/product/map/ 

2 http:/ /www. eso.org/science/healpix 



where < F >j corresponds to the mean emission com- 
puted for the set of pixels i, and < Nhi >i to the mean 
HI column density for the same set of pixels. Note that 
all the data-sets used here are cosecant-law subtracted 
(see Sect- IHS) . 

To keep high signal-to-noise ratio, only 5 sets of pixels 
are considered here, with increasing Nhi- The first set 
(labeled "0" in Eq. Q}), serves as the "reference" set 
and corresponds to the lowest column density regions 
(representing ~ 5% of the sky). We are thus left with 4 
sets of pixels k with increasing Nhi and derive accord- 
ingly four mean spectra F v (k). The sets of pixels are 
selected on the cosecant-law removed HI emission that 
can be negative. For reference, the total mean HI column 
density (i.e. non cosecant-law subtracted) for the 4 bins 
are 3.3, 4.1, 5.6 and 9.9 1 20 at/cm 2 . By construction, 
the spectra are normalised to 10 20 at/cm 2 . Note that F 
in Eq. represents alternatively the DIRBE, FIRAS, 
DMR, WMAP, free-free and synchrotron data. 
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Fig. 1. Spectrum of the Hl-correlated component normalised to 10 20 at/cm 2 (with increasing N(HI) from a to d). The 
black diamonds at 100 and 140 fim are the DIRBE data, the black curve is the FIRAS spectrum, the black crosses with 
error bars are the DMR data and the black triangles at 3.2, 4.9, 7.3 and 9.1 mm are the WMAP data. Also displayed 
are the free- free and synchrotron contributions (blue stars and green triangles respectively). The green continous line 
is the result of a fit of the DIRBE 100, 140 ^m and FIRAS spectra (200 < A < 500 /xm) with a modified Planck curve 
with a v 2 emissivity law (the so-called stable thermal dust component). The residual WMAP emission (which is the 
WMAP- free-free - synchrotron - stable thermal dust component) is shown as red stars. 




2000 



10000 2000 



10000 



Lambda (,um) 



Lambda (/im) 



Fig. 2. Zoom on the millimeter part of Fig. Symbols and colors are the same as on Fig. ^ Added is the residual 
WMAP emission after having removed only the free-free and synchrotron contributions (light-blue diamonds). 
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Component 


3.2 mm 


4.3 mm 


7.9 mm 


9.1 mm 


WMAP 


1 


2.05 10" 


12 


6.15 10" 


L3 


2.92 10" 


13 


2.19 10" 


13 




2 


1.52 10" 


12 


4.65 10" 


L3 


2.25 10" 


13 


1.80 10" 


13 




3 


7.59 10" 


L3 


2.26 10" 


L3 


1.31 10" 


13 


1.27 10" 


13 




4 


7.22 10" 


L3 


2.03 10" 


L3 


1.34 10" 


13 


1.37 10" 


13 


Free- free 


1 


1.46 10" 




1.02 10" 


13 


7.23 10" 


14 


6.03 10" 


14 




2 


1.14 10" 


L3 


7.93 10" 


L4 


5.65 10" 


14 


4.71 10" 


14 




3 


8.33 10" 


L4 


5.79 10" 


L4 


4.12 10" 


14 


3.44 10" 


14 




4 


7.36 10" 


L4 


5.12 10" 


L4 


3.64 10" 


14 


3.04 10" 


14 


Synchrotron 


1 


5.32 10" 


L!> 


4.28 10" 


L4 


6.56 10" 


14 


7.60 10" 


14 




2 


3.57 10" 


L4 


3.95 10" 


L4 


5.67 10" 


14 


7.24 10" 


14 




3 


4.27 10" 


L4 


3.40 10" 


L4 


5.01 10" 


14 


6.76 10" 


14 




4 


3.79 10" 


L4 


4.64 10" 


L4 


6.19 10" 


14 


8.20 10" 


14 


Stable thermal 


1 


3.17 10" 


L3 


3.95 10" 


14 


5.53 10" 


15 


1.86 10" 


15 


dust 


2 


2.81 10" 


L3 


3.50 10" 


L4 


4.90 10" 


15 


1.64 10" 


15 




3 


2.91 10" 


L3 


3.62 10" 


L4 


5.08 10" 


15 


1.71 10" 


15 




4 


3.66 10" 


L3 


4.55 10" 


L4 


6.39 10" 


15 


2.15 10" 


15 


Residue 


1 


1.58 10" 


12 


4.31 10" 


13 


1.48 10" 


13 


8.11 10" 


14 




2 


1.09 10" 


12 


3.12 10" 


L3 


1.02 10" 


13 


5.89 10" 


14 




3 


3.41 10" 


L3 


9.75 10" 


L4 


3.42 10" 


14 


2.29 10" 


14 




4 


2.45 10" 


L3 


5.94 10" 


L4 


2.95 10" 


14 


2.20 10" 


14 



Table 1. WMAP, free-free, synchrotron and stable thermal dust component emission (in W/m 2 /sr, normalised to 
10 20 at/cm 2 ) in the four HI bins together with the residual emission (which is equal to WMAP - Free-free - Synchrotron 
- Stable thermal dust). The mean total HI column densities are 3.3, 4.1, 5.6 and 9.9 10 20 at/cm 2 for the bin 1, 2, 3 
and 4 respectively. 



3.2. Results 

The four spectra are presented on Fig. with the DIRBE 
data points at 100 and 140 /Ltm, the FIRAS spectra 
(displayed only between 210 and 1000 yum), the DMR 
data points at 90, 53 and 31 GHz and the WMAP data 
points at 3.2, 4.9, 7.3 and 9.1 mm (all these data points 
are in black on Fig. 1). A zoom on the millimeter part 
of the figure is presented on Fig. |3 We fit the DIRBE 
100, 140 ^m and FIRAS spectra (200 < A < 500 fim) 
with a modified Planck curve with a v 2 emissivity law 
(the result of the fit is displayed on Fig. ^ and [2J). We 
know that this fit is inconsistent with FIRAS data below 
^500 GHz where an excess component is detected (Reach 
et al. 119951 Finkbeiner et al. Q2H9$- However, discussing 
this component is not the goal of this paper. We only 
concentrate on the millimeter part on the spectra and 
how it relates to the far-infrared emission. It is important 
to note that this far-infrared dust emission has a stable 
spectrum, not changing with increasing opacity (Lagache 
et al. 1999). In this framework, the v 2 modified black 
body is well representative and useful for comparison 
between spectral far-infrared and millimeter shapes 3 . 
This far-infrared dust emission extrapolated at millimeter 
wavelengths will be called the "stable thermal dust 
component" . 

First, we see on Fig.|5|that there is a strong millimeter 
excess (with both DMR and WMAP data) with respect 



to the stable thermal dust component (i.e the v 2 modified 
black body). This excess decreases significantly (by a 
factor of about 5 at 3.2 mm) when the HI column density 
increases, although the far-infrared emission remains 
nearly constant (at the ~6% level). The far-infrared 
emission is dominated by the so-called Large Grain 
dust component. The millimeter excess, which changes 
rapidely with opacity, is thus not likely associated with 
this dust component. 

We can go further by removing to the WMAP emission 
the corresponding free-free, synchrotron and stable ther- 
mal dust component contribution. The residual WMAP 
emission is shown on Fig.n an dEl(red stars) and detailed 
in Tabled First, at each frequency, the residual emission 
exhibits a strong decrease (by about a factor of 5) with HI 
column densities (from bin 1 to 4). Second, the residual 
emission decreases from 3.2 to 9.1 mm in each HI bin. On 
Fig. |21 are shown the WMAP residual emissions for the 4 
bins at 3.2, 4.9, 7.3 and 9.1 mm, normalised to the 90 GHz 
DMR residual emission (the 31 and 53 GHz DMR residual 
emissions have also been computed but are not displayed 
to avoid confusion. Results, although more noisy, are in 
very good agreement with WMAP). This figure shows that 
we do not detect any significant variations in the spectral 
shape of the residual emission 4 . Thus, the Hl-normalised 
residual emission, although decreasing in amplitude with 
the HI column density, has a constant spectrum. 



3 The way we are fiting the far-infrared stable component is 
not critical since we focus on the variable millimeter emission. 



4 This however will have to be quantified when smoothed 
WMAP data with error bars will be available. 
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HI bin 1 


HI bin 2 


HI bin 3 


HI bin 4 


uR v (3.2 mm) 


1.58 10~ 12 


1.09 10~ 12 


3.41 W^ 13 


2.45 10 -13 


uR v (4.9 mm) 


4.31 10 - " 


3.12 lO - " 


9.75 10~ i4 


5.94 10~ 14 


vR v (7.3 mm) 


1.48 lO"" 


1.02 lO - " 


3.42 10~ 14 


2.95 10~ 14 


uR v (9.1 mm) 


8.11 10~ 14 


5.89 10 -14 


2.29 10~ 14 


2.20 10~ 14 


i/I„(240 |m) 


1.81 10" s 


1.63 10" s 


1.59 10 _s 


1.83 10" s 


vI v (U0 fj,m) 


3.90 10~ s 


3.54 10 _s 


3.41 10~ s 


3.90 10 _s 


i/I„(100 pim) 


2.61 10" s 


2.54 10" s 


2.31 10" s 


2.43 10" s 


z;I„(60 (im) > 3° 
^I,vsg(60 pm) > 3° 
^I„(60 /im) \p\ > 15° 
^I„vsg(60 Mm) > 15° 


1.16 10" 8 
9.24 10" 9 
9.34 10~ 9 
7.01 10" 9 


9.37 10" 9 
6.98 10" 9 
7.23 10" 9 
4.83 10" 9 


7.18 10" 9 
4.94 10" 9 
4.57 10~ 9 
2.32 10~ 9 


6.48 10" M 
4.21 10" 9 
4.92 10" 9 
2.65 10" 9 



Table 2. Excess residual emission (z^R^) at 3.2, 4.9, 7.3 and 9.1 mm with 240, 140, 100 associated brightnesses (in 
W/m 2 /sr, normalised to 10 20 at/cm 2 ) for the 4 HI bins. The 60 Aim brightness is given for > 3° and |/3| > 15° to 
show that the decrease may not be due to any residual zodical light emission. Also given are the 60 Aim brightnesses 
corresponding to the Very Small Grains dust component only (we have removed from the 60 a™ emission, 1(60 Aim), 
the best v 2 modified black body fit done on the Large Grain dust component). 



4. Discussion 

To account for the galactic energy emitted from the 
mid-infrared to the submillimctcr, it is necessary to have 
a broad dust size distribution from large grains down to 
large molecules. For example, Desert et al. (1990) (see 
also Draine & Anderson 119851 Puget et al. 119951 Weiland 
et al. 119861 Siebenmorgen & Krugel 1992J Dwek et al. 
1997] and more recently Li & Draine l2"()()l(l have proposed 
a consistent interpretation of both the infrared emission 
in diffuse HI clouds and the interstellar extinction curve 
using a model with three components: PAHs (Policyclic 
Aromatic Hydrocarbons), Very Small Grains (VSGs) 
and Large Grains. PAHs and VSGs are small enough 
(a < 10 nm) to experience significant temperature 
fluctuations after photon absorption. They emit over a 
wide range of temperatures and dominate the emission 
for A < 60 /im. The Large Grain component is the more 
traditional dust component historically inferred from 
optical studies. These grains are in equilibrium with the 
incident radiation field with a temperature of about 17 K 
in the diffuse atomic medium (Boulanger et al. I1996|l . 
The Large Grain dust component is expected, at long 
wavelengths, to be proportional to the total amount of 
solid material. The large spatial variations of the infrared 
spectrum over the wavelenghts range 12-60 Aim have been 
interpreted as changes in the abundance of small grains 
(Boulanger et al. 119901 Laureijs et al. 119911 Bernard 
et al. 119931 Abergel et al. 11994(1 . In particular, strong- 
deficits of the transiently heated grains emitted at 60 Aim 
are observed in dense interstellar clouds, these deficits 
being explained by grain-grain coagulation processes (e.g. 
Stepnik et al. l2T70ljl . 

It has been shown in Sect 3.2 that the Hl-normaliscd 
residual WMAP emission (i.e. the excess above free-free 
and synchrotron contributions and the stable thermal 
dust component) is well traced at large scale by the 



3 p 



2 r 




12 3 4 

Bin number 

Fig. 3. Hl-normalised residual WMAP emission at 3.2 
(red), 4.9 (green), 7.3 (blue) and 9.1 (magenta) mm. All 
spectra are normalised on the DMR 90 GHz residual emis- 
sion for the first HI bin (black points with error bars) . The 
HI column density increases from bin 1 to 4 (from 3.3 to 
9.9 1 20 at/cm 2 ). No significant variations in the spectral 
shape of the residual emission is detected. 



HI gas and (1) exhibits a constant spectrum from 3.2 
to 9.1 mm but (2) significantly decreases in amplitude 
when Nhi increases, contrary to the far-infrared emission 
(associated with the so-called stable thermal dust com- 
ponent) which stays rather constant (cf. Table 2). It is 
thus very likely that the residual WMAP emission is not 
associated with the Large Grain dust component. The 
decrease in amplitude ressembles in fact to the decrease 
of the PAH/ VSGs emission observed in dense interstellar 
clouds. By extrapolating the PAH/ VSGs behaviour 
from dense interstellar clouds to the diffuse medium, we 
may expect, when increasing the HI column density, to 
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decrease the PAH and VSGs proportion and thus the 
mid-infrared emission. If this is true, then the PAH/VSGs 
Hl-correlated emission should decrease with HI column 
densities. This decrease, if present, is very hard to observe 
in the mid-infrared due to the strong residual interplane- 
tary dust emission at large scale. On the DIRBE bands, 
only the 60 /jm may be used. We have computed for the 
4 HI bins the 60 /xm Hl-correlated emission with two 
different cuts in ecliptic latitude (\/3\ > 3° and > 15°). 
Although the absolute level of the 60 /xm Hl-correlated 
emission varies for the 2 cuts, we observe nearly the 
same significant decrease of the 60 /mi emission with 
the HI column density (cf. Table 2). The 60 fim band 
may be contaminated by the Large Grain emission (30 
to 40 %, e.g. Desert et al. 1990). Therefore, we remove 
to the 60 /xm emission the Large Grain contamination 
using the best v 2 modified black body fit (Fig. 1). The 
observed decrease at 60 /xm becomes even larger (Table 
2). Although this result should be interpreted with care 
due to the zodiacal contamination at 60 /xm, it suggests 
that the WMAP residual emission is associated with the 
small transiently heated particles. 

The previous results suggest the anomalous microwave 
component is associated with the transiently heated dust 
particles, but its exact physical mechanism remains to be 
found. On the possible models of the anomalous emission 
linked to the transiently heated particles are: 

— The "spinning-dust" which is the rotational emission 
from very small dust grains (Draine & Lazarian, 1998a 
1998b). However, although the spinning dut emission is 
in good agreement with the WMAP emission at 7.3 and 
9.1 mm, it is inconsistent with the 3.2 mm emission. 

— The VSGs long-wavelength emission. VSGs are tran- 
siently heated when an ultraviolet photon is absorbed. 
The mean interval between successive ultraviolet pho- 
tons is longer than the cooling time and thus, between 
2 impacts, the temperature of the particles is very low 
(but is at least the CMB temperature). Such particles 
could therefore emit significant emission in the mil- 
limeter channels. 

The models have large uncertainties linked to the un- 
known properties of the small particles. It is therefore very 
difficult to predict the exact contribution of the two in the 
millimeter. 
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